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ABSTRACT 

We investigate the relative sensitivities of several tests for deviations from Gaussianity in the 
primordial distribution of density perturbations. We consider models for non-Gaussianity that 
mimic that which comes from inflation as well as that which comes from topological defects. 
The tests we consider involve the cosmic microwave background (CMB), large-scale structure 
(LSS), high-redshift galaxies, and the abundances and properties of clusters. We find that 
the CMB is superior at finding non-Gaussianity in the primordial gravitational potential (as 
inflation would produce), while observations of high-redshift galaxies are much better suited 
to find non-Gaussianity that resembles that expected from topological defects. We derive a 
simple expression that relates the abundance of high-redshift objects in non-Gaussian models 
to the primordial skewness. 

Key words: cosmology: theory - galaxies - clusters of galaxies - large scale structures - 
cosmic microwave background - methods: analytical 
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1 INTRODUCTION 

Now that cosmic-microwave-background (CMB) experiments (de 
Bernardis et al. 2000; Jaffe et al. 2000; Balbi et al 2000; Lange et al. 
2000) have verified the inflationary predictions of a flat Universe 
and structure formation from primordial adiabatic perturbations, 
we are compelled to test further the predictions of the simplest 
single-scalar-field slow-roll inflation models and to look for possi- 
ble deviations. Measurements of the distribution of primordial den- 
sity perturbation afford such tests. If the primordial perturbations 
are due entirely to quantum fluctuations in the scalar field respon- 
sible for inflation (the "inflaton"), then their distribution should be 
very close to Gaussian (e.g., Guth & Pi 1982; Starobinski 1982; 
Bardeen, Steinhardt & Turner 1983; Falk, Rangarajan & Srednicki 
1993; Gangui et al. 1994; Gangui 1994; Wang & Kamionkowski 
2000; Gangui & Martin 2000). However, multiple-scalar-field mod- 
els of inflation allow for the possibility that a small fraction of 
primordial perturbations are produced by quantum fluctuations in 
a second scalar field. If so, the distribution of these perturbations 
could be non-Gaussian (e.g., Allen, Grinstein & Wise 1987; Kof- 
man & Pogosyan 1988; Salopek, Bond & Bardeen 1989; Linde & 
Mukhanov 1997; Peebles 1999a; Peebles 1999b; Salopek 1999). 
Moreover, it is still possible that some component of primordial 
perturbations are due to topological defects or some other exotic 
causal mechanism (Bouchet et al. 2000), and if so, their distri- 
bution should be non-Gaussian (e.g., Vilenkin 1985; Vachaspati 
1986; Hill, Schramm & Fry 1989; Turok 1989; Albrecht & Steb- 
bins 1992). Detection of any non-Gaussianity would thus be invalu- 
able for appreciating the nature of the ultra-high-energy physics 



that gave rise to primordial perturbations. Ruling such exotic pos- 
sibilities in or out will also be necessary to test the assumptions that 
underly the new era of precision cosmology. 

There are several observables that can be used to look for pri- 
mordial non-Gaussianity. CMB maps probe cosmological fluctu- 
ations when they were closest to their primordial form, and many 
authors have developed various mathematical tools to test the Gaus- 
sian hypothesis. The statistics of present-day large scale structure 
(LSS) in the Universe can also be used (e.g.. Coles et al. 1993; 
Luo & Schramm 1993; Lokas et al. 1995; Chodorowski & Bouchet 
1996; Stirling & Peacock 1996; Durrer et al. 2000; Verde & Heav- 
ens 2000). The properties and abundances of the most massive 
and/or highest-redshift objects in the Universe also contain pre- 
cious information about the nature of the initial conditions (e.g., 
Chiu, Ostriker & Strauss 1998; Robinson, Gawiser & Silk 1999; 
Robinson, Gawiser & Silk 2000; Willick 2000; Matarrese, Verde & 
Jimenez 2000 (MVJ); Verde et al. 2000b). In Verde et al. (2000a; 
VWHKOO), the relative sensitivities of the CMB and LSS to several 
broad classes of primordial non-Gaussianity were compared, and it 
was found that forthcoming CMB maps can provide more sensitive 
probes of primordial non-Gaussianity than galaxy surveys. Here 
we extend the results of that paper to include comparisons to the 
abundances of high-redshift galaxies as well as the abundance and 
properties of clusters. One of our original aims was to determine 
whether any of these probes would be able to detect the miniscule 
deviations from Gaussianity that arise from quantum fluctuations in 
the inflaton; unfortunately, we have been unable to find any. Never- 
theless, some detectable deviations from Gaussianity are conceiv- 
able with multiple-field models of inflation and/or some secondary 
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contribution to primordial perturbations from topological defects. 
We will follow VWHKOO and parameterize the primordial non- 
Gaussianity with a parameter that can be dialed from zero (cor- 
responding to the Gaussian case) for two different classes of non- 
Gaussianity. We will then compare the smallest value for the pa- 
rameter that can be detected with each of the different approaches. 



2 THE METHOD 

2.1 Models for primordial non-Ganssianity 

There are infinite types of possible deviations from Gaussianity, 
and it is unthinkable to address them all. However, we can con- 
sider plausible physical mechanisms that produce small deviations 
from the Gaussian behavior and thus analyze the following two 
models for the primordial non-Gaussianity (e.g.. Coles & Barrow 
1987, VWHKOO, MVJ). In the first model, we suppose that the frac- 
tional density perturbation 5{x) is a non-Gaussian random field that 
can be written in terms of a Gaussian random field 4>{x) through 
(Model A) 

5 = + eA(<^'-(0')). (1) 

In the second model, we suppose that the primordial gravitational 
potential 3>(x) is a non-Gaussian random field that can be written 
in terms of a Gaussian random field ^(a) through (Model B) 

$ = + eB(</.' -(</.')). (2) 

Non-Gaussianity in the density field is then obtained from that in 
the potential through the Poisson equation. Here, $ and 5 refer to 
the primordial gravitational potential and density perturbation, be- 
fore the action of the transfer function that takes place near matter- 
radiation equality. 

Although not fully general, these models may be considered as 
the lowest-order terms in Taylor expansions of more general fields, 
and are thus quite general for small deviations from Gaussianity. 
The scale-dependence of the non-Gaussianity in the two models 
differs. Model A produces deviations from Gaussianity that are 
roughly scale-independent on large scales, while Model B produces 
deviations from non-Gaussianity that become larger at larger dis- 
tance scales. Although we choose these models essentially in an ad 
hoc way, the non-Gaussianity of Model B is precisely that arising 
in standard slow-roll inflation and in non-standard (e.g., multifield) 
inflation (Luo 1994; Falk, Rangarajan & Srednicki 1993; Gangui 
et al. 1994; Fan & Bardeen 1992; see also below). Model A more 
closely resembles the non-Gaussianity that would be expected from 
topological defects (e.g., VWHKOO). In either case, the lowest- 
order deviations from non-Gaussianity (and those expected generi- 
cally to be the most easily observed) are the three-point correlation 
function (including the skewness, its zero-lag value) or equivalently 
the bispectrum, its Fourier-space counterpart. It is straightforward 
to calculate these quantities for both Models A and B. 



2.2 Cosmic Microwave Background and Lat^e Scale 
Structure 

Temperature fluctuations in the CMB come from density perturba- 
tions at the surface of last scatter, so the distribution of temperature 
fluctuations reflects that in the primordial density field. It is thus 
straightforward to relate the density-field bispectra of Models A 
and B to the bispectrum of the CMB. Density perturbations in the 
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Figure 1. Afmax as a function of redshift. At a given redshift one should 
only consider those masses (< Mmax) for which at least one object is 
expected in the whole sky for Gaussian initial conditions. The shaded re- 
gion encloses predictions for Mmax(z) from different mass functions in 
the literature; we adopted the currently favored cosmological model with 
parameters: Slo = 0.3, Ao = 0.7, h = 0.65, crs = 0.99 and transfer 
function of Sugiyama (1995) with Q.^ = O.OlS/ft^ (ACDM). 

Universe today grew via gravitational infall from primordial per- 
turbations in the early Universe, and this process alters the mass 
distribution in a calculable way. Cosmological perturbation theory 
allows the bispectrum for the mass distribution in the Universe to- 
day to be related to that for the primordial distribution. 

VWHKOO calculated the smallest values of eA and ee that 
would be accessible with the CMB and with LSS. For the CMB cal- 
culation, it was assumed that a temperature map could be measured 
to the cosmic-variance limit only for multipole moments £ < 100; 
it was assumed (quite conservatively) that no information would be 
obtained from larger multipole moments. The LSS calculation were 
made under the very optimistic assumption that the distribution of 
mass could be determined precisely from the galaxy distribution 
(i.e., that there was no biasing) in a survey of the size of SDSS 
and/or 2dF. VWHKOO found that the smallest values of e that can 
be detected with the CMB under these assumptions is ea ~ 10~^ 
and eb ~ 20 (Komatsu & Spergel 2000 including noise and fore- 
ground but neglecting dust contamination found that es > 5 from 
the Planck experiment ), while the smallest values measurable with 
LSS are eA ~ 10"^ and eb ~ 10^. More realistically, the galaxy 
distribution will be biased relative to the mass distribution, and this 
will degrade the sensitivities to nonzero ea and eb obtainable with 
LSS. VWHKOO thus concluded that the CMB will provide a keener 
probe of primordial non-Gaussianity for the class of models consid- 
ered. 

2.3 High-redshift and/or massive objects 

According to the Press-Schechter theory, the abundance of high- 
redshift and/or massive objects is determined by the form of the 
high-density tail of the primordial density distribution function. 
A probability distribution function (PDF) that produces a larger 
number of > 3a peaks than a Gaussian distribution will lead to 
a larger abundance of rare high-redshift and/or massive objects. 
Since small deviations from Gaussianity have deep impact on those 
statistics that probe the tail of the distribution (e.g. ?, MVJ), rare 
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high-redshift and/or massive objects should be powerful probes of 
primordial non Gaussianity. The number densities of high-redshift 
galaxies and/or of clusters (at either low or high redshifts) provides 
a very sensitive probe of the PDF. Since the Gaussian tail is decay- 
ing exponentially at higher densities, even a small deviation from 
Gaussianity can lead to huge enhancements in the number densi- 
ties. 

The non-Gaussianity parameters ea.b are effectively "tail en- 
hancement" parameters (c.f., MVJ)^ 

In order to determine the minimum value of ea.b that can 
be detected using high-redshift objects, one needs to compute by 
how much the observed number density of objects changes with re- 
spect to the Gaussian case when the primordial field is described 
by equations |^ and (^). We calculate this enhancement using the 
results for the mass function for mildly non-Gaussian initial con- 
ditions obtained analytically in MVJ. Conservatively, we make the 
assumption that objects form at the same redshift at which they are 
observed (zc — 2); since for some objects the dark halo will have 
collapsed before we observe them, the assumption gives a lower 
limit to the amount of non-Gaussianity. 

The directly observed quantity, however, is not the mass func- 
tion, but is A'^(> M, z), the total number of objects -in the survey 
area- of mass > M that collapse at redshift z. In fact it is extremely 
difficult to obtain an accurate estimate of the mass of high-redshift 
objects, what is a more robust quantity is the minimum mass that 
these objects must have in order to be detected at that redshift. This 
quantity is related to the mass function, n{M, z), by 

/>oo 

N{>M,z)= / n{M,z)dM . (3) 

J M 

In calculating the enhancement of high-redshift objects due 
to primordial non-Gaussianity, we restrict ourselves to consider, at 
any given redshift, only those masses, M < Afniax(2:) for which 
at least one object is expected in the whole sky for Gaussian initial 
conditions (N{> A/max, z) = I'm An radians This is illustrated 
in Fig. 1 for a ACDM model (hereafter we adopt the currently fa- 



* In fact, when looking on a particular scale, it is always possible to param- 
eterize the deviation of the PDF from Gaussianity, with some "effective" eA 
or . if the PDF is not too non-Gaussian. It is easy to understand this state- 
ment if one thinks in terms of skewness. Physical mechanisms that pro- 
duce non-Gaussianity generically produce non-zero skewness in the PDF 
for the simple reason that underdense region cannot be more empty than 
voids while overdense regions can become arbitraiily overdense. Skewness 
can be scale dependent, but for a given value of the skewness there is one- 
to-one correspondence to ca.b parameters (see the Appendix), 
t This choice for the threshold N{'> Mmax,^) = 1 is motivated by 
the following considerations. Of course it is not robust to detect a non- 
Gaussianity that suppresses the number of objects with respect to the Gaus- 
sian prediction, since one can always argue that one did not look hard 
enough, or that the objects are there but are somewhat "invisible". So we 
set to detect a non-Gaussianity that enhances the number of objects rela- 
tive to the Gaussian case. If within Gaussian initial conditions we expect 
N{> M, 2) ~ in the whole sky, and observations find N{> M, z) > 1 
in the survey area, we can say that we have detected non-Gaussianity. How- 
ever, the non-Gaussianity (or tail enhancement) parameter is directly re- 
lated to the ratio of observed Nng{> M,z) to the Gaussian predicted 
N{> M, z) (see Eq. 4). Obviously this ratio is well defined for any 
A^ng > and N > 0, but the observed Nng can only be an integer > 1. 
The tail enhancement parameter will then make Nng > N (and we con- 
sider only cases where Nng > ION). It is reasonable therefore to consider 
only those masses and redshift for which the theoretical prediction for the 
Gaussian N is> 1. 



vored cosmological model with parameters: Q,o — 0.3, Aq — 0.7, 
h — 0.65, (Tg = 0.99 and transfer function of Sugiyama (1995) 
with fli, = 0.015//i^) where the shaded region encloses predic- 
tions for i\fnmx(2) from different mass functions (e.g.. Press & 
Schechter 1974; Lee & Shandarin 1998; Sheth & Tormen 1999; 
Jenkins et al. 2000). 

Given the rapidly dying tail of the Gaussian PDF, small un- 
certainties in the mass determination of high-redshift objects could 
lead to overestimate the value of ea.b. An overabundance of galax- 
ies of estimated mass Afe, which in principle can be attributed to a 
non-zero value of eA,B, can also be explained under the hypothesis 
of Gaussian initial conditions if the actual galaxy mass A/truc is 
A/truo < Me- We thus include conservative values for the uncer- 
tainty AM in the mass determination of high-redshift objects and 
we then calculate the minimum change AA'^ in the number density 
of objects over the Gaussian case that cannot be attributed to the un- 
certainty in the mass determination. For a given uncertainty in the 
mass, this can be computed by using the standard Press-Schechter 
(PS) theory (Press & Schechter 1974). Observationally it is difficult 
to measure the mass of high-redshift clusters with accuracy better 
than 30%, with either weak lensing or the X-ray temperature, and 
of high-redshift galaxies better than a factor 2 (AM = M; at least 
of their stellar mass). Although the calculations in this section are 
obtained using the standard PS theory, our conclusions will be es- 
sentially unchanged if we had used modified PS theories (e.g., Lee 
& Shandarin 1998; Sheth & Tormen 1999; Sheth, Mo & Tormen 
1999; Jenkins et al. 2000; see below). 

With the mass uncertainties discussed above, we obtain that 
the minimum AA'^ that cannot be attributed to AM is a factor 10 for 
clusters and a factor 100 for galaxies (see, e.g.. Fig. 6 of MVJOO). 

We therefore estimate the minimum eA.B that can be measured 
from the abundance of high-redshift objects as the one that corre- 
sponds respectively to a factor 100 and 10 change in the observed 
number density of objects (N{> M, z)) over the Gaussian case. 
This condition can be written as 



N„g{> M,z)/N{> M,z) = R{M,z) > R,, 



(4) 



where A' is obtained using the Gaussian mass function while Nng 
is obtained using the non-Gaussian mass function as in MVJ, and 
i?, is set to be 100 for galaxies and 10 for clusters. 

For small primordial non-Gaussianity (i.e., for small values 
for ea.b), it is possible to derive an expression for R{M, z) using 
the analytical approximation for the mass function n„g found in 
MVJ. Doing so we find 



R(M, z) 
Here, 



J^{aMM)-^eM-%f^\nM,z.,tK,^)dM 

M 

/.7(-MM)-iexp[-%i]|^|dA^ 



.(5) 



F{M, 2e,eA,B) 



dS: 



3,M 



6-\/l — S3,m3c{zc)/3 



dM 



— S3^mSc{Zc) /3 d(JM 



(JM 



dM 



and 



5t{Zc) = 5c(Zc)\/l - Ss,mSc{Zc)/3, 
Sc{Zc) = Ac/D{Zc), 

where D{zc) is the linear-theory growth factor, and Ac is the linear 
extrapolation of the over-density for spherical collapse. 



(6) 

(7) 
(8) 
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In the formulas above, Ss^m denotes the primordial skewness, 



c (1) / ^ 

03,A/ = eA.BMs.Af/'^M 



(9) 



where the expressions for ^3 and aj^j can be found in MVJ sec- 
tion 3.2 equations (37) and (38). However, for Sz > l/5c(zc), 
the mass function n„g (M, z) has to be evaluated numerically and 
equation is not valid. 

For the cosmological model considered here and the redshifts 
of interest, the quantity Ac takes a nearly constant value (~ 1.686) 
in the PS theory. A better fit to the mass function of halos in high- 
resolution N-body simulations is however obtained by lowering Ac 
for rare objects and giving it an extra mass and redshift dependence 
(Sheth & Tormen 1999; Bode et al. 2000), as motivated by ellip- 
soidal collapse (e.g., Lee & Shandarin 1998; Sheth, Mo & Tormen 
1999). 

It is possible to understand the effect of a lower Ac by the 
following argument. For rare fluctuations such as high-redshift ob- 
jects one is probing the mass function above the knee. Since the 
mass function drops very rapidly as M increases we can approx- 
imate A'^(> M, Zc) ~ n(M, Zc)M. It is then possible to obtain 
an analytic expression for r(M, 2c) = nng{M, Zc)/n{M, z^) ~ 
R{M, Zc) if the primordial non-Gaussianity is small: 



r(M, 



;exp 



S3 3, 
3 



'dcrjv 



1 - 



SaSc 



.(10) 



For a given mass M, r{M, z^) slowly decreases when low- 
ering Ac, slightly damping the effect of non-Gaussianity. For ex- 
ample when lowering Ac from the value we assume here 1.686, to 
the value « 1.5 — appropriate to fit the numerical mass function of 
Sheth & Tormen (1999) for the range of masses and redshifts con- 
sidered here — r(M, Zc) decreases by less than a factor 2. However 
this effect is compensated by the fact that, by lowering Ac, ob- 
jects are created more easily also with Gaussian initial conditions, 
and it is therefore possible to consider objects of higher M and/or 
z, where the effect of non-Gaussianity is bigger. In summary, the 
conclusions obtained by assuming Ac = 1.686 will not be sub- 
stantially modified. 

It is important to note that for Model A, the primordial skew- 
ness has the same sign as ea, while for Model B the primordial 
skewness has the opposite sign of that of eb. In detecting non-zero 
eA,B from CMB maps, the sign of the skewness does not influ- 
ence the accuracy of the detection of non-Gaussianity, but, when 
using the abundance of high-redshift objects the sign of the skew- 
ness matters. Only a positively skewed primordial distribution will 
generate more high-redshift objects than predicted in the Gaussian 
case. Although a negatively skewed probability distribution will 
generate fewer objects than the Gaussian case, a decrement might 
be difficult to attribute exclusively to a negatively skewed distribu- 
tion. Therefore in the following we will consider only negative eb 
and positive ea- 



2.3.1 Cluster size-temperature distribution 

Verde et al. (2000b) showed that the size-temperature (ST) distri- 
bution of clusters is fairly sensitive to the degree of primordial non- 
Gaussianity. If clusters are created from rare Gaussian peaks, the 
spread in formation redshift should be small and so should the scat- 
ter in the ST distribution. Conversely, if the probability distribution 
function has long non-Gaussian tails, then clusters of a given mass 



we observe today should have a broader formation redshift distri- 
bution and thus a broader ST relation. In Verde et al. (2000b), the 
non-Gaussianity considered is a log-normal distribution; it is not 
strictly equivalent to Models A or B considered here. However, for 
small deviations from Gaussianity, the two models can be identi- 
fied if, for a given scale, they produce the same skewness in the 
density fluctuation field. We thus find that in the ACDM model the 
minimum ea and eb detectable with the ST distribution method are 
3x10^"^ and 500 respectively. These estimates assume that the cos- 
mology and (78 are well known, but use only the local cluster data 
set of Mohr et al. (2000). Of course, with improved observational 
data, the ST method could probably yield stronger constraints. 



3 RESULTS 

We find that the non-Gaussianity of Model A has a bigger effect 
on high-redshift galaxies than on high-redshift clusters. This can 
be understood for the following reason. For Model A the skewness 
5*3, A/ is approximately scale independent (dSa^M /dM = 0). Thus, 
as found in MVJ, the mass function for non-Gaussian initial con- 
ditions is obtained from the PS mass function for Gaussian initial 
conditions replacing Sc{zc) — > S,{zc). The effect of a non-zero 
skewness is therefore to lower the effective threshold for collapse 
thus allowing more objects to be created. For a given 5*3, <5*(2c) is 
a monotonically decreasing function of Zc. Since galaxies can be 
observed at Zc much bigger than that of clusters, the effect is big- 
ger. On the other hand, clusters are better probes than galaxies for 
Model B. In fact, for Model B the induced skewness in the den- 
sity field is scale dependent and the effect of non-Gaussianity is 
roughly the same for galaxies with 8 < 2 < 10 and clusters with 
1 < 2 < 3. However since mass determinations are more accurate 
for clusters than for galaxies, we have i?,, dusters < -R», galaxies: 
clusters are therefore better probes. 

In Fig. I we show the ratio R = iV„<,(> M, z)/N{> M, 2) 
(cf. eq. (5)) for galaxies at redshift 2 = 8, 9 and 10 for ea = 
5 X 10^* (model A, left panel) and clusters at redshift 2 = 1,2 and 
3 for eb = —200 (model B, right panel), as a function of M. Lines 
are plotted only for masses where, for Gaussian initial conditions, 
one would expect to observe at least one object in the whole sky 
with the most conservative estimate (see Fig. |l]). Note that those 
high-redshift objects represent 3- to 5-a peaks. If we now require 
R{M, Zc) > 7?,, we deduce that the minimum detectable deviation 
from Gaussian initial conditions will be ea 5 x 10"'* (from high- 
redshift galaxies) and |eB | ~ 200 (from high-redshift clusters). We 
also estimate that an uncertainty of 10% on as would propagate 
into an uncertainty of 25% in eb (from clusters) and of 70% in ea 
(from galaxies). 

The minimum es detectable from high-redshift cluster abun- 
dances is much larger than the value that can be measured from the 
CMB (eb ~ 5 to 20 for Planck data), while for ea, high-redshift 
galaxies are much better probes than the CMB, which can only de- 
tect eA ~ 10"^. 

We therefore conclude that ;/ future NGST or 30- to 100-m 
ground-based telescope observations of high-redshift galaxies yield 
a significant number of galaxies at 2 ~ 10 and are able to deter- 
mine their masses within a factor 2, these observations will perform 
better than CMB maps in constraining primordial non-Gaussianity 
of the form of Model A with positive ea- Conversely, forthcoming 
CMB maps will constrain deviations from Gaussianity in the ini- 
tial conditions much better than observations of high-redshift ob- 
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Figure 2. Ratio R{M, z) = N„g{> M, z)/N{> M, z) for galaxies at redsliift z = 8,9 and 10 for ea = 5 X 10""' (left panel) and clusters at redshift 
z = 1, 2 and 3 (right panel), for eg = 200 as a function of M. Lines are plotted only for masses where, for Gaussian initial conditions, one would expect to 
observe at least one object in the whole sky with the most conservative estimate (see Fig. Note that these high-redshift objects represent 3- to 5-cr peaks. 
The values for the number density enhancement R that can safely be attributed to primordial non-Gaussianity are R = 100 for galaxies (left panel) and 
i? = 10 for clusters (right panel). See text for details. 



jects for Model B (with positive and negative value for eb) and for 
Model A with negative ea- 



3.1 Slow-roll parameters and primordial skewness 

The type of non-Gaussianity of Model B is particularly interesting 
because initial conditions set from standard inflation show devia- 
tions from Gaussianity of this kind. In fact, it is possible to relate 
the two slow roll parameters, 



IGtt 



Tj* = 



8n 



V 



(11) 



to the non-Gaussianity parameter eb- In equation ( |l l| ) rnpi is the 
Planck mass, V denotes the inflaton potential and V' and V" the 
first and second derivatives with respect to the scalar field. The 
skewness Ss for $b, 53,$ = ($b)/{'1?b)^' can be evaluated fol- 
lowing a similar calculation of Buchalter & Kamionkowski (1999), 
obtaining 



53,* = 2eB X 3[1 + 7(r- 



(12) 



where 7(n) <C I and weakly depends on n if n < 0, but diverges 
for n > Q. For a scale-invariant matter density power spectrum, 
n — —3, 7(n) = 0, and so 5*3,* = 6eB- 

We can then compare this expression with the value for the 
skewness parameter for the gravitational potential arising from in- 
flation to infer the magnitude of eb- Gangui et al. (1994) calculate 
the CMB skewness for the Sach- Wolfe effect S2 in several infla- 
tionary models; 52 is related to Sa,* by S2 ~ Ss^is.A'^^ where 
Asw = 1/3. From this it follows that S2 ~ 353,$ — ISeb- The 
condition for slow roll from Gangui et al. 1994 is ^2 < 20; thus, 
eb < 1, and the relation to the slow-roll parameters is (cf., Wang 
& Kamionkowski 2000) 



EB = (5/2)6, - (5/3)r?. 



(13) 



Since this combination of the slow-roll parameters is different from 
the combination that gives the spectral slope n of the primordial 
power spectrum (n — 2e, — 6r;, + 1), in principle, if eb could 
be measured with an error 1, it would be possible to determine 



observable 


min. |eA 


min. |eB 


CMB 


10-3 ~ 10-2 


20 


LSS 


10-2 


10^ 10* 


High z obj . 


(+)5 X 10-4 (gal.) 


(-) 200 (clusters) 


ST relation 


(+)3 X 10-3 


(-) 500 



Table 1. Minimum | e a | and | eb I detectable form different observables and 
their sign when positive skewness is required for detection. For Model A the 
primordial skewness has the same sign as ep^, while for Model B the pri- 
mordial skewness has the opposite sign as eg- In detecting non-zero eA,B 
from CMB maps, the sign of the skewness does not influence the accu- 
racy of the detection of non-Gaussianity, but, when using the abundance of 
high-redshift objects it is robust to detect non-Gaussianity that produces an 
excess rather than a defect in the number density. Only a positively skewed 
primordial distribution will generate more high-redshift objects than pre- 
dicted in the Gaussian case. 



the shape of the inflaton potential through eq. (jllj). However, from 
the present analysis, an error of eb of about an order of magnitude 
larger seems to be realistically achievable. 



4 DISCUSSION AND CONCLUSIONS 

We considered two models for small primordial non-Gaussianity, 
one in which the primordial density perturbation contains a term 
that is the square of a Gaussian field (Model A), and one in which 
the primordial gravitational potential perturbation contains a term 
proportional to the square of a Gaussian (Model B). The non- 
Gaussianity of Model B is precisely that arising in standard slow- 
roll inflation and in non-standard inflation, while Model A more 
closely resembles the non-Gaussianity that would be expected from 
topological defects. We investigated the relative sensitivities of sev- 
eral observables for testing for deviations from Gaussianity: CMB, 
LSS and high-redshift and/or massive objects (e.g., galaxies and 
clusters). 

The analytic tools developed above allow us to address the 
question of whether the abundance of currently known high- 
redshift objects can be accommodated within the framework of in- 
flationary models for a given cosmology. Recently Willick (2000) 
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has studied in detail the mass determination of the cluster MS 1054- 
03 concluding that its mass lies in the range 1.4±0.3xl0^^ Mq for 
f2m = 0.3 (similar to the independent mass estimates by, e.g., Tran 
et al. (2000) and Newmann & Arnaud (2000)). As already pointed 
out by Willick (2000), for S7m > 0.3 the expected number of ob- 
jects like MS 1054-03 in the survey area is < 0.01; i.e., it must be 
a 3-(j fluctuation or larger. Using the formalism we have described 
here, a primordial non-Gaussianity parameterized by ee > 400 
would be required to account for MS 1054-03 as a Icr fluctuation in 
the ACDM model described above. This value is much too large to 
be consistent with slow-roll inflation. Our calculation shows that if 
such a non-Gaussianity exists, it would be easily detectable from 
forthcoming CMB maps. 
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APPENDIX 

In this Appendix we quote the expressions for the primordial bis- 
pectrum and skewness for the two non-Gaussian models considered 
in this paper. The large-scale structure (LSS) bispectrum for model 
A is (e.g. VWHKOO) 



B(ki, ka, ks) = 2eAP(fci)P(fc2) + eye. 



(14) 



where P denotes the power spectrum. The cosmic microwave back- 
ground (CMB) bispectrum for model A is (e.g. VWHKOO) 



B. 



(112(3 



(2^1 + l)(2i;2 + l) (2^3 + 1) ,t,i^t,. 

\o j 



4:TT 



, 2£A 

9 



2^ ^ 



(15) 



where Ce denotes the cosmic microwave background power spec- 
trum, g denotes the radiation transfer function and (...) denotes 
the Wigner 3J symbol. The LSS bispectrum for model B is (e.g. 
VWHKOO) 



B(ki,k2,k3 



P(fci)P(fc2)2eB 



MkiMk2 



+ eye. 



(16) 



where^A^fc ~ (2fc2r(fc)(l + .z))/(3//c?) and T denotes the matter 
transfer function. The CMB bispectrum for model B is (e.g. Luo 



t This expression is strictly valid only for an Einstein de Sitter Universe, 
for a more general model A4 is defined by 5^ (2) = Mk{^)^{k) where "I" 
denotes the gravitational potential field. 
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1994; Wang & Kamionkowski 2000,VWHK00,Komatsu & Spergel 
2000): 



(2£i + l) (2^2 + 1) (24 + 1) n,i,i. 



47r 



X — (C«iCf2 + eye 



(17) 



The corresponding primordial skewness 83 = {5"^) / {S'^)'^ 
where S denotes Sp/pfor the large-scale structure case and AT/T 
for the cosmic microwave background is easily obtained from the 
consideration that (5^) is given by: 



(27r)3 (27r)3 
(in the absence of spatial filtering) and 



d"fc3B(ki,k2,k3)cS"(ki+k2+k3)(18) 



- y 



(2^i + l)(2£2 + l)(2£3 + l) 
47r 



(19) 



for LSS and CMB respectively. For example in the large scale 
structure case, model A, for a power law power spectrum and in the 

absence of spatial filtering^ we have that Sz = 6eA- 



§ The expression for {&''^)lss 'he general case can easily be derived 
following the calculations of Buchalter & Kamionkowski (1999) by setting 
61 = and 62/2 = EA- 
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